function retval = Dm_coupling (input_Bo,input_H)
 
    global b d p g Rf OA         %%% 常数         
    global Bo II u ii            %%% 时间向量，方位角向量     
    global H W h Ds Dm Os Om     %%% 参数向量          
    global PointSet_three_Bo_u   %%% 三相线点集 
    global PointSet_down_Bo_u_O  %%% 气液面下部点集    
    global PointSet_up_Bo_u_O    %%% 气液面上部点集        
    %%% 特殊 全局变量
    global Point_r_three_Bo_u  Point_z_three_Bo_u  a_down_Bo_u       
    
    for N=1:II
     if input_Bo(1,N)~=0
         retval(1,N)=1.279*input_Bo(1,N).^(-0.176)*(input_H(1,N));
     else
         retval(1,N)=0;
     end
end
